Theory of optical transitions in graphene nanoribbons 
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Matrix elements of electron-light interactions for armchair and zigzag graphene nanoribbons are 
constructed analytically using a tight-binding model. The changes in wavenumber (An) and pseu- 
dospin are the necessary elements if we are to understand the optical selection rule. It is shown that 
incident light with a specific polarization and energy, induces an indirect transition (An = ±1), 
which results in a characteristic peak in the absorption spectra. Such a peak provides evidence that 
the electron standing wave is formed by multiple refiections at both edges of a ribbon. It is also 
suggested that the absorption of low-energy light is sensitive to the position of the Fermi energy, 
direction of light polarization, and irregularities in the edge. The effect of depolarization on the 
absorption peak is briefly discussed. 



I. INTRODUCTION 

The dynamics of the electrons in graphene is governed 
by an equation that is similar to the relativistic equation 
for the massless Dirac fermion»ii^ Since the Fermi veloc- 
ity of a carrier is about lO^m/s, the Dirac fermion moves 
forward 1 nm in 1 fs. Therefore, in nanometer-sized 
graphene, the carrier reaches the edge before its motion 
is affected by perturbations such as electron-phonon and 
electron-electron interactions. As a result, the electronic 
properties of the system are sensitive to the presence of 
an edge. The behavior of electrons near the edge of a 
graphene sample is unique due to the reflection of the 
massless Dirac fermion. 

In graphene nanoribbons, the importance of the edge is 
marked by reflections of electron taking place at the both 
edges of the ribboui^"— The reflections result in the forma- 
tion of a standing wave of a Dirac fermion. In this paper, 
we examine optical transitions in graphene nanoribbons 
and clarify characteristic features of the standing wave of 
a Dirac fermion. Knowing the rules for the optical tran- 
sitions is an important step in understanding the optical 
properties of graphene nanoribbons. 

Graphene edges are categorized into two groups; arm- 
chair and zigzag edges with respect to the symmetry 
of the hexagonal lattice^"— It is known that standing 
waves near the armchair edge and near the zigzag edge 
are distinct for various reasons, such as their pseudospin 
and Berry's phase.l^ In this paper, we study nanorib- 
bons with armchair and zigzag edges in great detail, and 
briefly discuss the effect of irregularities in the edge on 
the absorption spectra. 

Here, we mention previously published literature on 
the optical absorption of graphene nanoribbons. Hsu and 
Reichiii investigated the absorption of linearly polarized 
light parallel to a zigzag nanoribbon and found that a 
direct transition is not allowed, in contrast to the case 
of nanotubesJ^ Gundra and Shuklapi^ pointed out that 
the polarization dependence on the absorption spectra is 
important for characterizing the edge structure. Whereas 
these studies were based on numerical simulations, our 
study provides analytical results for electron-light matrix 



elements. The analytical result clearly shows the effects 
of pseudospin and momentum conservation on the optical 
properties of a graphene nanoribbon. 

This paper is organized as follows. In Sec. [Hi we deal 
with armchair nanoribbons. By constructing the matrix 
elements of the electron-light interaction, we show that 
dynamical conductivity depends on the direction of the 
polarization of the incident light with respect to the ori- 
entation of the edge. In Sec. IIIIl we examine absorption 
spectra for zigzag nanoribbons. The effect of edge irreg- 
ularities on dynamical conductivity is studied in Sec. IIVI 
Our discussion and conclusion are provided in Sees. IVl 
and IVIl respectively. 



II. ARMCHAIR NANORIBBON 

In this section we study the electron-light interaction 
in armchair nanoribbons. In Sec. Ill A] we review the elec- 
tronic properties of armchair nanoribbons to provide the 
necessary background. In Sec. IIIBl the matrix elements 
of electron-light interaction are constructed. The results 
are used to clarify the absorption spectra in Sec. Ill CI 



A. Electron wavefunction 

The energy dispersion relation for armchair nanorib- 
bonsii^iS is written as 



e^{k,0) = S7ov^l + 4cos2 6' + 4cos6'cosfc/, (1) 

where 70 is the hopping integral (70 = 3 eV) and I = 
V^a/2 (a is a lattice constant [a — 2.46 A]). The energy 
is characterized by the band index s and two parameters 
fc, and 9. The superscript, s, represents the conduction 
(valence) energy band and takes values s — 1 (—1) on the 
right-hand side of Eq. ([T]), fc is the wave vector parallel 
to the edge, and 6 stands for the phase in the direction 
perpendicular to the edge [see Fig.[T][a)]. By setting fc = 
and = 2tt/3 into Eq. (P), we see that £^^(0, 27r/3) = is 
satisfied. Thus, the conduction and valence energy bands 



touch at the point {k, 9) = (0, 2tt/3), which is called the 
Dirac point. The energy dispersion relation for armchair 
nanoribbons Eq. ([T|) is the same as that for graphene."i^ 
There are two independent Dirac points (known as K and 
K' points) in the Brillouin zone (BZ) of graphene, on the 
other hand, there is a single Dirac point in the BZ of 
armchair nanoribbons. That is, although e^{k,9) is zero 
at the other point {k,d) = (0,— 27r/3), this point is not 
included in the BZ of armchair nanoribbons. In fact, the 
BZ of armchair nanoribbons is given by G (0,7r) and 
kb e [— 7r,7r), where b = \/3a, as shown in Fig. [ijb). The 
BZ of armchair nanoribbons covers only one-half of the 
graphene's BZ because the reflection taking place at the 
armchair edge identifies 6 with 
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FIG. 1: (a) The structure of an A'' armchair nanoribbon. The 
unit length along the edge is denoted by b (= \^a). Carbon 
atoms are divided into A (•) and B (o) atoms, (b) The BZ of 
armchair nanoribbons: 6 £ (0, tt) and kb £ [— tt, tt). The circle 
stands for the Dirac point, (c) The definition of the polar 
angle 0{k, 6). The circle represents an energy contour (1 eV). 
(d) Schematic diagrams for intra- and inter-band transitions 
from the initial state with 6„ = (on the 6'-axis). 



The wave function is given by 
1 



(j)'j{k,6) = ^e-'''^iJ~i) sin J6 



-ie(kfi) 



(2) 



where J (= 1, . . . , N) is the coordinate perpendicular to 
the edge [see Fig. [lja)]i^"— A detailed derivation of the 
wave function is given in Ref. Il6l The upper (lower) com- 
ponent of Eq. ([2]) represents the amplitude at A-atom 
(B-atom) in the box shown in Fig. [TJa). The relative 
phase between the two components, 6(fc, 9) in Eq. ([2]), is 
the polar angle defined with respect to the Dirac point 
as shown in Fig. [TJc). The wave function with 9 = 27r/3 
does not have an amplitude at J = 3, 6, . . . due to sin( 76*) 



of Eq. ([2]) . The corresponding line nodes have been ob- 
served in recent scanning tunneling microscopy topogra- 
phy,— which evidences the wave function Eq. ([2|). 

The wave function vanishes at J = and J = N + 1 
(i.e., at fictitious edge sites). The boundary condition for 
J = is given by <pj^Q{k,9) — 0, which is satisfied with 
arbitrary values of 9. The phase 9 is quantized by the 
boundary condition for J = A'' -I- 1 , ^^_|_j(fc, 9) = 0, as 



N + 1' 



{n^l,...,N), 



(3) 



where n represents the subband index. Meanwhile, we 
assume that the wave vector k parallel to the edge is a 
continuous variable. Because k is not changed by the per- 
turbations, which preserve translational symmetry along 
the edge, hereafter, we abbreviate (j)j{k,9n), Q{k,9n), 
and e'^{k,9n) by omitting k and 9 as 
respectively. 



n 7) ©n, and £*, 



B. Selection rule 

The electron-light interaction is written as -ffcm = 
—ev ■ A, where — e is the electron charge, v — (v^jVy) 
is the velocity operator, and A ex e~''^*e is a spatially 
uniform vector potential. Here, tku corresponds to the en- 
ergy of the incident light and e denotes the polarization. 
The matrix elements of v are classified into inter-band 
{(t>n\v\(j)m) and intra-band {(|)'f^\v\(j)^) transitions. The 
contributions of these transitions to the optical property 
of a nanoribbon depend not only on huj and e, but also 
on the position of the Fermi energy s-p. The intra-band 
transition may be omitted only when ep = (band cen- 
ter) at zero temperature. In general, it is necessary to 
consider both transitions. 

First, we consider the intra-band transition. The ma- 
trix elements of v are calculated using Eq. ([2]). The de- 
tails of the calculation are provided in Appendix [X] The 
results are 

m — n E even 

{^'n^\v.\<f>'n) = { .2 Vp ^ ^^ ^ 

-I {crx)mn m-ne odd, 



n m — n 

{<Pm\Vy\(f>n) = SmnVF{cry)Tnn, 



(5) 



where 6mn is the Kronecker delta, v-p (= ^qI/K) is the 
Fermi velocity, and a^ (z = x, y) denotes the Pauli matri- 
ces. The matrix elements of cr, are written as 



lfe'e„^g-.e„ 



^yjynn — „ (6 



iBr, 



-ie„ 



(6) 
(7) 



(Ti is the pseudospin, which brings certain features to the 
optical transitions. For example, in Eq. ([6]), {ax)mn — 1 
is satisfied for forward scattering from 0„ = to Qm = 0, 
while {<Jx)mn = for backward scattering from 0„ = to 
©m — T^ [see Fig. [TJd)]. Since {ax)mn is proportional to 



((/)^|fa;|</>^) as described in Eq. (jj]), the x-polarized light 
(A^) does not cause intra-band backward scattering. 

In Eq. ((4]), the factor in front of the pseudospin {ax)mn 
arises from momentum conservation. It can be shown 
that momentum conservation for {4>^\vx\4>n) leads to the 
foUowing summation with respect to the out-of-phase 
trigonometric functions (Appendix E)) : 



2 
— ^ sin( J6l„i) cos{ J9n) - 

J=i 






in - 


- 71 e even 


2 1 


m - 


- n e odd, 


TT m — n 




(8) 



where 9m(n) is given by Eq. ([3]). The summation takes a 
non-zero value only when m — n {= An) is an odd num- 
ber. The momentum conservation for (0mkyl<^n) leads 
to the summation of the in-phase trigonometric functions 
[see Eq. (jA10|) ]. so that the summation takes a non-zero 
value only for m ~ n. Hereafter, we call the process sat- 
isfying m — n Cz odd an indirect transition for which the 
wavenumber of the initial state changes (An S odd), and 
the process satisfying m = n is called a direct transition 
(An = 0). 

We have seen that the velocity matrix element is deter- 
mined by two factors: pseudospin and momentum con- 
servation. Next, we consider the inter-band transition 
based on this understanding. The calculated matrix ele- 
ments are given by 







K) 



Vp 



TT m — n 



^ y I mn 



- n G even 

- n e odd. 



(9) 
(10) 



In Eq. (ITU)) . 5mn shows that the y-polarized light (Ay) 
results in a direct inter-band transition [see Fig. [Ijd)]. 
Thus, the transition amplitude depends on the diagonal 
matrix element of the pseudospin {ax)nn- From Eq. ([6]), 
we see that {ax)nn takes a maximum value of {ax)nn = 
cos 6„ = ±1 for 6„ = or tt. On the other hand, {ax)nn 
vanishes for &„ — ±7r/2. Therefore, the electrons on the 
6'-axis (8„ = or tt) are selectively excited by Ay, while 
the electrons near the fc-axis (6„ = ±7r/2) are excited 
very little. 

As we can see in Eq. ([9]), the x-polarized light {Ax) 
results in an indirect transition [see Fig. [T](d)]. An inter- 
band transition amplitude from 0„ = to Qrn = tt is 
enhanced because the strength of the pseudospin takes 
a maximum value of \{(Jy)mn\ = 1- This indirect inter- 
band transition is a forward scattering, namely, the elec- 
tron crosses the Dirac point. An inter-band transition 
that is not across the Dirac point, such as a (backward) 
transition from 0„ = to Qm = 0, is allowed by mo- 
mentum conservation if to — n is an odd number. How- 
ever, it is strongly suppressed by the pseudospin because 
{o'y)mn — for the process. 



C. Dynamical conductivity 

With the optical selection rule established by Eqs. (|4|), 
([S]), dH), and (fTOl) . let us investigate the absorption of 
light. The dynamical conductivity is given by (a = x,y) 

2 



h 






s' s nrn/c 



/[4J-/[4] (0^|(-ez;„)|0^) 






+ huj + i6) 



(11) 

where f[e] is the Fermi distribution function, S is the 
nanoribbon area, and S is inversely proportional to 
the relaxation time of the excited electrouiiS. The real 
part of the dynamical conductivity, Ke{axx) (R-e(crj,y)), 
represents the absorption spectrum of x-polarized {y- 
polarized) light. In Eq. pT|) . we assume ^ = 10 meV 
and room temperature when evaluating the Fermi distri- 
bution function. 

The solid curve in Fig. [5Ja) shows Re(o'j,j,) calculated 
for an TV = 15 armchair nanoribbon. The width of this 
nanoribbon is 1.7 nm. As we have shown in Eq. (jlOp . 
the absorption of y-polarized light results from a direct 
inter-band transition. Thus, the peaks denoted by A, 
B, and C of Re((Tj,j,) originate from the direct inter-band 
transitions < ^ </)^ (A), <^5;_i -^ (j)l_^ (B), and <+i ^ 
•^n+i (C), respectively, as shown in Fig.I^Jb). The energy 
position of the A peak corresponds to the energy gap 
(-£'gap) of 0.66 eV for an iV = 15 armchair nanoribbon. 
Both a large density of states at the band edge (fc — 0) 
of each subband and the pseudospin are important as 
regards enhancing these peak intensities. 

The dashed curve in Fig. [2Ka) represents Re{axx)- As 
we have seen in Eq. ©, the absorption of a:-polarized 
light is the result of an indirect inter-band transition. 
Thus, the peak denoted by D of Re{axx) is due to the 
transitions (p^ — )■ (p^-i and 4>n-i ~^ 4'n^ ^^ shown in 
Fig- IHc). Because the resonance energy (e^ — e^_i) is 



between e^ 



eZ and e'i 



which are the res- 



onance energies for the direct inter-band transitions A 
and B, the peak D appears between the A and B peaks. 
Note that another indirect transition 0^^ — >■ 0^_|_i is also 
allowed by momentum conservation. However, the pseu- 
dospin strongly suppresses the backward transition for 
the states at the band edge (fc = 0) and the large density 
of states does not result in an absorption peak. Hence, 
there is no prominent peak for Re{(Jxx) between the B 
and C peaks. Moreover, the peak corresponding to the 
lowest-energy backward transition is the most prominent, 
and the peak for the second lowest (and higher) energy 
transition, such as (/>^ — )■ (/>5;_3 and <J3^j^i — )■ 4>n-2^ is gen- 
erally suppressed by the momentum conservation giving 
rise to the suppression factor of (m — n)~^ in Eq. ([Q]). 
To clarify this point, we show the calculated Re(o'a;a;) for 
an iV = 55 armchair nanoribbon with i?gap — 0.2 eV 
in Fig. [3l[a). For Ke{axx), there is only a peak on the 
low-energy side. 

It is clear that the absorption spectrum of an arm- 
chair nanoribbon exhibits strong anisotropy for the po- 
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FIG. 2: (a) The calculated Re(craa) for an A'^ = 15 (semi- 
conducting) armchair nanoribbon. The inter-band transition 
processes relevant to the absorption peaks of Ke{ayy) and 
^e{axx) are shown in (b) and (c), respectively, (d) The opti- 
cally allowed indirect inter-band transition (from (jj^ to (j>'n-i 
or from 4>^-i to (f)"^). Bottom: The pseudospin (denoted by 
arrows) is useful for showing that only forward scattering is 
allowed in the inter-band optical transitions. 



larization direction of incident light. ^■^ More specifically, a 
number of peaks appear when the polarization is parallel 
to the armchair nanoribbon, while only a single promi- 
nent peak appears on the low-energy side when the po- 
larization is set perpendicular to the ribbon. Because 
the energy gap of (semiconducting) armchair nanorib- 
bons behaves as £^gap oc N~^^-~ the energy position of 
the peak for perpendicular polarization is inversely pro- 
portional to the width of the armchair nanoribbon. 

Since ep is generally not zero, it is meaningful to point 
out that the absorption of low-energy photons is sensi- 
tive to the position of s-p. To explain this feature, we 
show the dynamical conductivity for an iV = 55 arm- 
chair nanoribbon with ep = 0.3 eV in Fig. [Sjb). For 
-E-gap < ?Jw < 2|£p|, Re((Tj,y) is suppressed by the Pauli 
exclusion principle, while a peak remains for I{,e{axx)- 
This peak is attributed to forward intra-band scattering. 
In addition, the Drude peak at ?iuj = appears only for y- 
polarized light. This feature arises from the fact that the 
diagonal matrix element of Vy can take a non-zero value, 
while that of Vx is zero, as we have seen in Eqs. (U) and 



FIG. 3: Re(craa) of an A'^ = 55 armchair nanoribbon with an 
energy gap about 0.2 eV, calculated for (a) ep = eV and 
(b) £p = 0.3 eV. 
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III. ZIGZAG NANORIBBON 



In this section, we focus on zigzag nanoribbons. We 
show that the optical properties of zigzag nanorib- 
bons are quantitatively different from those of armchair 
nanoribbons. For example, when the polarization of an 
incident light is parallel to the zigzag nanoribbons, the al- 
lowed inter-band transitions are indirect transitions. On 
the other hand, when the light polarization is perpen- 
dicular to the ribbons, a direct inter-band transition is 
allowed. The behavior of zigzag nanoribbons contrasts 
with that of armchair nanoribbons. 



A. Electron wavefunction 

The energy dispersion relation for zigzag nanoribbons 
is given by 



£"(fc,6') =s7o^ 



4 cos^ ( — I 
2; 



4 cos ( — I COS0, 



(12) 

where k is the wave vector afong the zigzag edge, and 9 
denotes the phase in the direction perpendicular to the 
edge [see Fig. gl^a)]. By putting {k,9) = (47r/3a,0) and 
(27r/3a,7r) into Eq. ([12]), we have e*(fc,6') = 0. The two 
Dirac points are non-equivalent, i.e., the BZ of zigzag 
nanoribbons contains the two independent Dirac points 
(K and K') [see Fig. UJ^b)]. This contrasts with the fact 
that the BZ of armchair nanoribbons has a single Dirac 
point. The difference in the number of Dirac points is 
because the reflection of an electronic wave at the zigzag 
edge constitutes intravalley scattering, while the reflec- 
tion at the armchair edge is intervalley scattering. 
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FIG. 4: (a) The structure of an N zigzag nanoribbon. (b) 
The BZ of a zigzag nanoribbon: 6 G (0, tt) and ka £ [0, 27r). 
The empty circles show the positions of the K and K' points. 
The dashed semicircles around the K and K' points represent 
an energy contour (1 eV). (c) The plot of 9n{k) for an A'^ 
zigzag nanoribbon. 



The wave function is expressed by 



^yik,9) = ^le^-'<^(' 



2VN 



-ie(k,e) 



-ije / e 



,j:e(fc,e) 



(13) 

where J (= 1,...,N) is a coordinate perpendicular to 
the zigzag edge [see Fig. Hlja)]. A detailed derivation 



of the wave function including the pseudospin is given 
in Ref. |22|. The wave function Eq. (|13p is a standing 
wave formed by the superposition of two waves propagat- 
ing in opposite directions. The first term on the right- 
hand side of Eq. (|13p represents an incident wave, and 
the second term corresponds to an edge reflected wave. 



The phase e(fc, 9) is defined through \e%k, 6l)|e-*®('='^) = 
2cos(fca/2) -I- e~'^r^ by which we can see that Q{k,9) 
corresponds to the polar angle defined with respect to 
the K or K' point. For the K point, Q{k,9) is given as 
shown in Fig.Ul^c). Note that the signs in front of Q{k, 9) 
for the constituent waves are opposite. This feature can 
be understood in terms of the reflection of pseudospin, 
namely, the zigzag edge alters the direction of the pseu- 
dospin of an incident wave.-- By contrast, the armchair 
edge does not change the direction of the pseudospin. In 
fact, in this case, we can reproduce the wave function for 
armchair edge Eq. ^ as 



JJd / e 



,-ie 



-i.i0 / e 



,-ie 



oc sm{J9) 



-40 



(14) 



It is the characteristic feature of the zigzag edge that 
the reflection of the momentum of the incident wave 
{9 -^ —9) accompanies the change of the pseudospin 
{Q — >■ —Q). This correlation between momentum and 
pseudospin does not exist for armchair nanoribbons. 
The wave function has two components 



^J 



Va.j 



s[n{J9-e{k,9)) 
s sm{ J9) 



(15) 



where the flrst (second) component represents the ampli- 
tude at A-atom (B-atom) [see Fig. Hl^a)]. Note that the 
minus sign in front of the second term on the right-hand 
side of Eq. (|13p ensures that the wave function satisfles 
the boundary condition at J = given by ip^ j^q = 0. 
The phase 9 is quantized as 9n by the boundary con- 
dition at J = N + 1, (f^ .j=N+i = 0- This boundary 
condition gives the constraint condition between k and 
9n, 6(fc, 9n) — {N + l)9n = —nTT, which is rewritten as 






(16) 



We plot 9n with n = 0, 1 as a function of k in Fig. SJc), 
where 0„ is a curved rather than a straight lineJ^ This 
feature also contrasts with that of an armchair nanorib- 
bon, for which the quantized 9n is independent of k [see 
Eq. ([3])]. Note that the curved line of 0„(fc) with n — 
reaches the fc-axis near the Dirac point. It can be shown 
that 9o{k) acquires an imaginary part after crossing the 
/c-axis as 9o{k) = i/£^{k). In Eq. (TTSt . £,{k) corresponds to 
the localization length. By putting 9o{k) — i/£,{k) into 
Eq. pBj) . we see that Qo{k, 9o{k)) also acquires an imag- 
inary part as Qo{k,9o{k)) = i{N + l)/^{k). This means 
that the pseudospin is polarized by the localization. The 
localized states near the zigzag edge are known as the 
edge stateSf2,i2, The curved line is essential to the existence 
of the edge localized states in the zigzag nanoribbons. 



B. Selection rule 

We start by showing the velocity matrix elements with 
respect to the inter-band transition (see Appendix [B] for 
derivation) , 



(<Pmbal<P'n) 






(17) 



The calculation for the pseudospin of the standing wave, 
on the right-hand side of Eq. p7)) , is straightforward with 
the use of Eqs. (fT3|) and (|16p . The results are 



{'^liWxWn) =-sine 



(7t+m)7r-|-(en+e„) 
2 



sm 



(7i-m)7r+(e„-e„,) 
2 



{n + m)TT + (6„ + 6,„) {n - ■m)n + (6„ - 6,„) 
cos Qm f sin [{n + m)7r -|- (0„ + 9m)] sin [{n — m)7r + (8„ — Qj. 



{n + m)TT + (e„ + 9™) 



{n - ■m)n + (9„ - 9,„) 



+ (nOm), (18) 



Hy^ml'^vl'Pn) 



sin9r 



sm 



(K-|-m)7r+(e„+e„) 
2 



sm 



-t- 



(n~m)7r+(e„-e„.) 
2 



(n -t- m)7r -f (9„ + 9m) (n - m)7r -f (9„ - 9™) 
cos 9,„ f sin [{n + mjn + (9„ 4- 9„,.)] sin [{n — m)-n + (9„ — 9,„)] 



(n + m)7r + (9„ + 9™) 



(n — m)7r + (9^ 



9. 



— (n o to). (19) 



r 



First, we consider a case where the polarization of the 
incident light is parallel to the zigzag nanoribbon. By 
setting TO = n in Eq. ((TO)) , we have ((y9^Jcry|(/?^) — 0. Since 
((^$^1(7^1^^) = leads to {•^%,\v,A^l) = in Eq. m. we 
conclude that the x-polarized light does not cause a direct 
inter-band transition. Therefore, the possible inter-band 
transition is an indirect one. To further explore the inter- 
band transitions produced by the x-polarized light, let us 
examine the electrons with 9 — ti/2 (or k = 47r/3a). By 
putting 9„ — tt/2 and 9m = t^/2 into Eq. (|T9)) . we obtain 



i^r. 



\Vn 



VF- 



2 sin^ (f Ar 



(20) 



Equation ^U^ suggests that the a;-polarized light induces 
indirect transitions when m~n (= An) is an odd number. 
The inter-band transitions with An = ±1 have advantage 
over the transitions with An ~ ±3, 5, . . . in producing 
prominent peaks in the dynamical conductivity since the 
suppression due to the momentum conservation is mini- 
mum when An = ±1. 

Next, we consider a case where the polarization of the 
incident light is perpendicular to the zigzag nanoribbon. 
When m — n, Eq. (fTSl) leads to 



{VnWx\(p'r, 



COS 9„ 



sin 9„ 

rnr + 9, 



(21) 



The second term on the right-hand side is suppressed by 
the factor of {TTn)~^ when n > 1. For n = 0, {(p'f^\ax\<fn) 
is zero when 9o = and —1 when 9o = tt. Thus, the y- 
polarized light gives rise to direct inter-band transitions. 



by which the electrons near the fc-axis (9 = or tt) are 
selectively excited. Since there is not a large density of 
states for the states near the fc-axis, the direct inter-band 
transition does not result in a prominent absorption peak. 



By putting 9,; 
obtain 



7r/2 and 9m = 7r/2 into Eq. (HH]), we 



/ c I I m ■ 2sin^ 



[f(n- 



1)] 



(n 



1) 



(22) 



Equation (1221) suggests that the ^/-polarized light gives 
rise to an indirect inter-band transition when n -I- m is 
an even number. For example, the transition ipQ -^ ip^ 
[which corresponds to the peak C in Fig. [5fc)] satisfies 
this condition. However, the amplitude of this process is 
small due to the suppression by momentum conservation. 
The indirect inter-band transitions caused by y-polarized 
light can be neglected. 

We have shown that the a;-polarized light (parallel to 
the zigzag edge) results in indirect inter-band transi- 
tions (An = ±1) for the states near the band edges.— 
The j/-polarized light induces direct inter-band transi- 
tions (An — 0) selectively for the states near the fc-axis. 
The polarization dependence of the inter-band optical 
transition in zigzag nanoribbons exhibits a 90° phase 
shift with respect to that derived for armchair nanorib- 
bons. 
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FIG. 5: (a) Calculated Re(crQQ) for an A'' = 11 zigzag 
nanoribbon. The optical transition processes relevant to (b) 
Re{axx) and (c) Re((Tyy). The flat band at the band center 
represents the energy dispersion of the edge states. 



The velocity matrix elements for the intra-band tran- 
sition are written as 



From Eqs. ^^ and ([IT]), we see that 



(23) 



(24) 



The equations mean that the selection rule for intra-band 
transitions can be obtained from that for inter-band tran- 
sitions by changing the polarization direction. For the 
intra-band transitions, the a;-polarized light results in a 
direct transition (An — 0), while the y-polarized light 
results in an indirect transition (An = ±1). 



C. Dynamical conductivity 

In Fig. mja), we show the calculated Re{aaa) for an 
A^ = 11 zigzag nanoribbon. The width of the ribbon is 
2.2 nm. The solid curve represents Ke{axx), where the 
A and B peaks are attributed to the transitions shown 



in Fig. ini^b). The change in the wavenumber, which is 
relevant to the A peak, is An = ±1, while that relevant 
to the B peak is An — ±3. The intensity of the B peak 
is much smaller than that of the A peak due to sup- 
pression by momentum conservation. The dashed curve 
represents Re(crj,j,), where the small C peak is due to the 
transitions shown in Fig. [SJc). The finite conductivity 
at low energy below the C peak originates from the di- 
rect transition shown in Fig. [Sjc). The conductivity is 
suppressed near huj = 0. This feature can be explained 
by the pseudospin of the edge states. Later we discuss 
the fact that the velocity matrix element between two 
edge states vanishes because the edge states are pseu- 
dospin polarized states (i.e., eigen state of az), so that 
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FIG. 6: The calculated Ke{aaa) of an N = 50 zigzag nanorib- 
bon for different Fermi energies: (a) £f = eV and (b) 
£F = 0.3 eV. 

Note that the indirect inter-band transitions A, B, and 
C in Fig. IHb) and (c) , involve the states near the Dirac 
point. Hence, we can expect the dynamical conductivity 
to be sensitive to the position of ep. To elucidate this 
point, we show the calculated Ke{(Taa) for an A^ = 50 
zigzag nanoribbon with different ep in Fig. [H^a) [ep — 
eV] and (b) [ep = 0.3 eV]. In Fig.IHKa), the most promi- 
nent peak appears for Ile{axx)j which corresponds to the 
transition A (</?q -^(fl) shown in Fig. EKb). The other 



peaks are due to (/jJ^ — )• (p^+i '^ith J^ = 1, 2, By con- 
trast, at a low energy in Fig. IH^b), the peak intensity 
is suppressed for Re(cra;x) and the most prominent peak 
appears for Tie{ayy). This change of the polarization de- 
pendence is because intra-band transitions, such as D 
shown in Fig. [5][c), are activated by y-polarized light. 
Finally, the Drude peak appears only when the polar- 
ization is set parallel to the zigzag nanoribbon, because 
the x-polarized light provides a direct intra-band transi- 
tion. This polarization dependence of the Drude peak is 
analogous to the case of armchair nanoribbon. 



IV. IRREGULAR EDGE 

In this section, we study the effect of edge irregularities 
on the dynamical conductivity by employing a numeri- 
cal simulation based on a tight-binding model. Although 
the variety of irregular edge structures that we consider 
is quite limited, the results help us to understand the way 
in which the absorption spectra are changed by irregu- 
larities. 

Figure[7{a) displays the calculated Re{(Taa) for a defec- 
tive A'^ = 15 armchair nanoribbon (thick lines). We intro- 
duced irregularities consisting of zigzag edges along the 
armchair edge of the ribbon. One of the zigzag parts is 
shown in the right-hand side of Fig. ^a). The size of the 
zigzag part is about Aa; = 5 A and Ay — 2 nm. For com- 
parison, we show Ke{aaa) for a pure armchair nanorib- 
bon without any irregularities (thin lines) in Fig. [7l[a). 
These results show that the irregularities do not induce 
any significant change in the ^e{axx) of the regular arm- 
chair nanoribbon. A notable effect of the irregularities 
appears for Re((Tyy) as the new peak at 0.4 eV, denoted 
by Z. Since the original pure armchair nanoribbon has 
an energy gap of 0.66 eV, the Z peak originates from ir- 
regularities consisting of the zigzag parts. Note that the 
appearance of the Z peak accompanies a reduction in the 
peak intensity at 0.66 eV. 

Because a small segment of the zigzag edge about 1 nm 
in length creates a zero-energy edge state in the energy 
spectrum,^ it is reasonable to consider that the Z peak 
is relevant to the edge state caused by the irregularity. 
Moreover, the reduction in the peak intensity at 0.66 eV 
suggests that the optical transition from the edge state 
to the state at the band edge of the first subband is ac- 
tivated. This consideration is consistent with the fact 
that the Z peak position (0.4 eV) is approximately half 
of the energy band gap of the pure armchair nanoribbon 
(0.66/2 eV). 

An interesting point here is that the edge state exist- 
ing near the zigzag irregularity can contribute to the op- 
tical transition, on the other hand, the edge state in pure 
zigzag nanoribbons does not. That is, the edge state near 
a zigzag irregularity is not identical to the edge state in 
a pure zigzag nanoribbon. For a pure zigzag nanoribbon, 
the edge atoms at J = are all A- atoms, while those 
at J = N are all B-atoms [see Fig. Hl^a)]. As a result. 
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FIG. 7: The calculated Ile{aaa) in the presence (thick lines) 
or absence (thin lines) of irregularities for (a) A'^ = 15 arm- 
chair and (c) A'' = 11 zigzag nanoribbons. (b) compares 
Re(o-aa) in the presence of an irregularity consisting of two 
sublattices (thin lines) with that consisting of one sublattice 
(thick lines) . In the simulation, we assume that the length of 
the nanoribbon is 40 nm, 5 = 50 meV, and ep = eV. 



the wave function of the edge state localized at J = is 
written as 



fACce M j ' 



(25) 



while the wave function of the edge state localized at 
J = N is written as 



_«^ /O^ 
V^B oc e « K . , 



(26) 



where ^ is the localization length. Note that these edge 
states have an amplitude only on the A-atoms or B-atoms 



depending on the edge atom. Since the electron-light ma- 
trix element is proportional to the matrix element of a^ 
or <7y, the matrix elements with respect to the pseudospin 
polarized edge states vanish as 



(1 0)a, 



= (i = x,j/). 



(27) 



One might consider that the matrix element of ai be- 
tween the edge state localized at J = and that localized 
at J = A^ can be non-zero because 



(0 



1 



^0. 



(28) 



However, such matrix elements are also strongly sup- 
pressed because ^ < iV is satisfied and there is little 
spatial overlap between the two edge states. Therefore, 
the edge state in a pure zigzag nanoribbon does not con- 
tribute to the optical transition. 

By contrast, the edge atoms of the zigzag irregularity 
include equal numbers of A-atoms and B-atoms, and the 
A- and B-edge atoms are close in space. Thus, the edge 
state at the A-edge atoms [(/3A(r)] can interfere with the 
other edge state at the B-atoms [(/3B(r)] to form a new 
edge state pair [(^^'^(r) = (pA(r) ± V&{y)\. This feature 
of the interference between two sublattices is essential to 
the optical transition. This consideration is consistent 
with the result shown in Fig.[71Jb), where the intensity of 
the Z peak decreases when the zigzag edge part consists 
only of one sublattice<^ Note also that the inter-band 
transition (p'"{v) — >■ (/9°(r) for the edge states contributes 
to the absorption spectrum at zero-energy [Tilj « 0) as 
shown in Fig. ^a). 

In Fig. IDJc), we show 'R,e{aaa) for a defective A^ = 11 
zigzag nanoribbon (thick lines) and that for a pure one 
(thin lines). The irregularities consist of armchair edges. 
One of the armchair parts is shown in the right-hand 
side of Fig. ^c). The irregularities do not result in a 
notable change of the absorption spectrum. This feature 
can be rationalized by using the fact that the armchair 
edge does not change the pseudospin via the reflection of 
the incident electron wave [see Eq. ([T4|]. 

The numerical results presented in this section suggest 
that for a defective armchair nanoribbon, the absorp- 
tion spectrum of transversely {x) polarized light is robust 
against the presence of irregularities, while the absorp- 
tion spectrum of low-energy longitudinally [y) polarized 
light is sensitive to them. For a defective zigzag nanorib- 
bon, the absorption spectrum appears insensitive to the 
irregularities. 



V. DISCUSSION 

The optical selection rule of armchair nanoribbons re- 
minds us of the optical selection rule of single- wall car- 
bon nanotubes (SWNTs),i? The selection rule of SWNTs 
shows that light polarized parallel (perpendicular) to the 



axis results in a direct (indirect) inter-band transition 
satisfying An = (An = ±1). It is known that the in- 
direct transition originates from the cylindrical topology 
of SWNTs. Although the indirect transition is interest- 
ing, the corresponding absorption peak is suppressed by 
the depolarization effectJ^ Indeed, the experimental re- 
sults for SWNTs do not show absorption peaks relevant 
to the indirect inter-band transition. This is consistent 
with the presence of a depolarization fieldi^"— While 
the topology of nanoribbons is quite different from that 
of SWNTs, there is a possibility that a depolarization 
field suppresses the absorption of the x-polarized light 
for armchair nanoribbons. The selection rule for zigzag 
nanoribbons is distinct from that for armchair nanorib- 
bons, and the zigzag nanoribbons seem to be anomalous 
with respect to the depolarization effect. The depolar- 
ization field points along the j/-axis, which is orthogonal 
to the direction of the polarization of the incident light. 

The optical transition between two subbands n and m 
depends on changes in the wavenumber An [— m — n) 
and the pseudospin {ai)mn- These factors An and {(Ji)mn 
could be treated independently in a discussion of the op- 
tical transitions for armchair nanoribbons. On the other 
hand, the factors mix and cannot be treated indepen- 
dently for zigzag nanoribbons. The pseudospin can be 
changed by interactions that are not included in our anal- 
ysis. For example, a potential difference between the A 
and B sublattices is induced by the stack in the Bernal 
configuration and gives rise to a polarization of pseu- 
dospin for the states near the Dirac pointJ^ Further- 
more, the pseudospin can be related to real spin because 
the edge states in zigzag nanoribbons are prone to be 
spin polarized if we take into account the Coulomb in- 
teractioni^i^ Although the selection rules that we have 
obtained provide the basis for exploring the optical phe- 
nomena exhibited by graphene nanoribbons, it is possi- 
ble that these interactions modify the absorption peaks 
present in this paper J^ 



VI. CONCLUSION 



We have paid considerable attention to the indirect 
optical transitions satisfying An — ±1. An inter-band 
transition with An = ±1 is induced when the polariza- 
tion of the incident light is perpendicular (parallel) to 
the armchair (zigzag) edge, while an intra-band indirect 
transition is induced when the light polarization is per- 
pendicular to both the armchair and zigzag edges. The 
indirect optical transitions can be traced to the inter-edge 
electronic coherence. The signal of the indirect transition 
appears as a low-energy peak in absorption spectra. The 
corresponding absorption peak remains in the presence 
of the irregularity and with different Fermi energy po- 
sitions, while it may be suppressed by a depolarization 
field. 



10 



Acknowledgments 



Eq. (|M]) . we have 



K.S. is supported by a Grant-in-Aid for Specially Pro- 
moted Research (Grant No. 23310083) from the Ministry 
of Education, Culture, Sports, Science and Technology. 



Appendix A: Derivations of Eqs. (Q, ((5]), ([9]), and 

GHD 

In this appendix, we give the derivations of the ma- 
trix elements of the velocity operator v for armchair 
nanoribbons. The operator v is defined as the first 
derivative of the nearest-neighbor tight-binding Hamil- 
tonian, H{A), with respect to a vector potential. A, as 
—ev = dH{A)/dA\A=o- The electron-light interaction 
is given by Hem — —ev ■ A, and the matrix elements are 
written as 

,7=1 

(,5/i+G+C„7+i + ^/i./C,7 + Sh-G'<f>lj_,) . (Al) 

On the right-hand side of Eq. (JA1[) . we have defined the 
translational operators. 



G+ = 



1 



G- ^ ( J ^_i, ) , (A2) 



and the 2x2 matrices of the potential induced by A as 

6hj = o-_A • Ri — cr+A • Ri, 

(5/i+ = <T_A-R2-cr+A-R3, (A3) 



5h 



J 



a-A ■ R3 - (j+A • R2, 



where 



0-+ 



1 
0/ ' 





1 



(A4) 



In Eq. (jAll) . the translational operators G^ should be 
inserted between 5h^ and 0^ jj^-^ because 0^ j of Eq. ^ 
is defined along the diagonal line of the boxes shown 
in Fig. [IJa). Although G^ ~ ctq for the states near 
the Dirac point, we do not use such an approximation 
here and keep the derivation mathematically exact. In 
Eq. (jA3|) . Rq (a = 1,2,3) are the vectors pointing from 
an A-atom to the nearest-neighbor B-atoms defined as 
follows: 



n/3 



R2 
R: 



1 



~Ltp(-;Go- ^ Qinc.Gli 



_ -^3 1 

3 c\ '-^cc^a; " TT'^cc^y: 



(A5) 



where e^; (e^) is the dimensionless unit vector for the x- 
axis (y-axis) and ace is the bond length between nearest- 
neighbor carbon atoms. By putting Eq. (jA5|) into 



5hj = -i(Ty{Ayacc), 

Sh+ ^ — —cr^iA^a^c) + -cTyiAyacc), (A6) 

Shj = —cTx{Axacc) + -cTyiAyacc). 

It is possible to consider a case where A is spatially mod- 
ulated.^^ But, here we assume a spatially uniform vector 
potential A, so that the potential does not depend on J. 
By putting Eqs. (PSj) and (KE\i into Eq. (IXTj) . and by 
comparing the coefficient of A^ (Ay) on both sides of the 
equation, we obtain 

TV 

N 



*™KI'/':)=«FV[C„7]M'^.-*f^, 



7=1 ^ "^ 



(A7) 
where the matrix T (T^^) is defined by T ee a^G^ 



T = ( i i ] , T-i 



Jkb 



ikb 



e- 

1 



(A8) 



Note that the eigen equation for the electron in an arm- 
chair nanoribbon is written in terms of these T, T^^ , and 
K^ = ax+ efjCTo asi2, 

TK,j+i + K^Ij + T- V:,7-i = 0. (A9) 

The energy eigen equation has been used to obtain 
(0f„|z;,|(/<^)ofEq. im. 

Below, we execute the summation over J in Eq. (|A7p . 
First, we consider ((?!>m|i'i,|0^) of Eq. (|A7p . By putting 
Eq. ([2]) into Eq. (IA7p . and by using the formula. 



we have 



2 ^ 
— ^ sin( 76*™) sin (J6'„) = 5,, 

7=1 



l^'yl'/'n) = VY5mn{c^y)mn, 



(AlO) 



(^m I ^y 



Cl^'ylO = il'Ff^mn <{ ((Ta;)m„ - y 



(All) 



These results are shown in Eqs. ^ and pH)) . In Eq. pUj) . 
the small term, eJi/3, has been omitted. Next, we take 
((/)^|wx|C) of Eq. dMJ. By using Eq. ^, we obtain the 



equation. 



(A12) 
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By putting this into Eq. (jA7|) . we obtain 



(d"xi0^) = ~WF 



2 sin 9„ 



.7=1 



y^ sin (J6lm) cos (J6>„) 
(A13) 



In this equation, we approximate 9n ~ 27r/3 or 
2sin0„/\/3 w 1 for the states near the Dirac point. 
Moreover, kl for those low-energy states is a small num- 
ber compared with the polar angles 6„ and 8„i. Thus, 
gg»(e™-fc/) _,. g-j(e„-fcO approximates to e*®" + e"*®", 



which is simply 2{ax)mn for s 
J in Eq. (|JT3t leads to 



2 ^ 
— 2_^ sm (J6',„) cos( J6'„) 

J=i 



1. The summation over 




771 — n G even 



+ 0{N-'^) m-neodd. 



(A14) 



The correction of ©(iV^^) to the matrix element has been 
neglected in the text. Thus, we reproduce Eqs. (jlj and 
([9]). To obtain the right-hand side of Eq. (jA14p . the quan- 
tization of Eq. ([3]) is essential. Namely, the presence of 
the armchair edge at J = iV plays an essential role. 



Putting Ra of Eq. (jA5l) into the above equations, we get 



7=2 (^9{k)A^ax + g{k)Ay<7y) a^c, 



Sh 

5h+ 

5h^j = —a-f-A- 



cr_ Ayflcc, 



y^cC' 



By inserting Eq. (jB3|) into Eq. (JB1|) . we have 



{'P'in\Vy\ip'„ 



(B3) 



J=i 

N . 2 1 

= WF ^ [(P™, j] ^ I -9{k)(Ty + i-e^a^ | (^^ j. 



(B4) 



where g{k) = 2cos(fca/2) and g{k) — (2/-\/3) sin (fca/2). 
Hence, g = —1 {g = 1) for the K (K') point, and 5 = 1 
for the K and K' points. To obtain the matrix elements 
of Vy , we utilized the energy eigen equation. 



Appendix B: Derivations of Eqs. (|17p and (|23p 

In this appendix, we give the derivation of the velocity 
matrix elements for zigzag nanoribbons. The electron- 
light matrix elements are given by 

N 

(^^J {-ev ■ A) 1^^) = -7ef Y. y^j] ^ X 



J=i 



(^'ijV,^i,7+i + Shjip^,j + Shjip',^j_^) 



where 



5hj = (7^A- e-*-R2 + e^-R3 



6h-, 



— a I A ■ 



cr A Ri, 



~R3 + e''~Ri 



-cr+A • Ri. 



(Bl) 



(B2) 



0'+V«„7-l +Kn(Pn,J- 



■^-<,7+i=0, (B5) 



where K^ = g{k)ax+e!^ and J = 1, . . . , A^. By multiply- 
ing az with the eigen equation from the left, we obtain 



^+<,7-i + <^^K<j - ^-flj+i = 0. 



(B6) 



Note that the sign in front of a^ becomes minus. This 
equation has been used in obtaining the matrix elements 
of Vy. In the text, the term proportional to e^ in Eq. (jB4p 
has been omitted. 
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